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Abstract 

Gaussian copulas are widely used in the industry to correlate two 
random variables when there is no prior knowledge about the co- 
dependence between them. The perturbed Gaussian copula approach 
allows introducing the skew information of both random variables into 
the co-dependence structure. The analytical expression of this copula 
is derived through an asymptotic expansion under the assumption of 
a common fast mean reverting stochastic volatility factor. This pa- 
per applies this new perturbed copula to the valuation of derivative 
products; in particular FX quanto options to a third currency. A 
calibration procedure to fit the skew of both underlying securities is 
presented. The action of the perturbed copula is interpreted compared 
to the Gaussian copula. A real worked example is carried out compar- 
ing both copulas and a local volatility model with constant correlation 
for varying maturities, correlations and skew configurations. 

1 Introduction 



Copula models arise in the market place when only quoted information about 
the behaviour of single assets is available but nothing or very little is known 
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about their joint relations. Products that depend on several assets allow to 
imply from the market some information about their joint behaviour. How- 
ever, there is a lot of market information embedded in just a single price and 
therefore, it is necessary to make some assumptions about the joint relations 
in order to imply some meaningful parameters out of market prices. On the 
other hand, when products which depend on several assets are not available 
in the market or they are not liquid enough, it is necessary to make some 
assumptions about the joint relations of their underlying assets in order to 
describe them in a simple and intuitive way through some parameters whose 
values might be given as input. Then, the dependence and sensitivity of 
a given product to these "unobserved" parameters allow to know the risk 
associated with them and allows to take a conservative position for trad- 
ing and managing them. Most of the assumptions of co-dependence among 
assets can be highly simplified through a copula model. One of the most 
popular contexts in which copula models have been used is credit (e.g. the 
valuation of Colateralized Debt Obligations). Another popular application 
appears in hybrid models which combine two different asset classes for which 
co-dependence information is not available. The particular application ad- 
dressed in this research work is FX (foreign exchange) quanto options to a 
third currency different from the two currencies of the underlying FX pair of 
the option. 

A copula model allows to obtain a joint probability distribution of two 
random variables when only their marginal distributions are known. As it is 
clearly explained in [5], the joint distribution of two random variables can be 
cleanly decomposed into two separate contributions: the distribution of the 
co-dependence (the copula function) and the marginal distribution of each 
random variable. This means that the copula function embeds completely 
the co-dependence information. The gaussian copula is the most popular, 
well-known and widely used copula model in the industry. It has become a 
reference proxy in the industry for it is analytical, easy to use and provides 
intuition. The main assumption of the gaussian copula is to consider that 
the random variables are normally distributed and their joint distribution 
is multi-variate normal. The copula function would only depend on the 
correlation matrix of the multivariate-normal distribution. This allows to 
embed the whole information of the co-dependence in just the correlation 
parameter. This hypothesis might not be reasonable when the distribution 
of the underlyings is skewed (the slope of the volatility surface with respect 
to strike is different from zero) or in the credit context, the probability of 
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extreme events (the tail of the gaussian distribution) is too low. 

The perturbed gaussian copula comes into play to improve the gaussian 
copula under the assumption that the two random variables have a common 
fast mean reverting stochastic volatility factor. This hypothesis allows to 
introduce the skew effect in the co-dependence of the two random variables. 
The joint distribution of the random variables is approximated through an 
asymptotic expansion to first order calculated using perturbation theory un- 
der the asumption that the common stochastic volatility factor is fast mean- 
reverting. This allows to obtain an analytical expression for the joint and 
marginal distributions of both random variables and therefore the copula 
funtion (the joint density divided by the product of the two marginal den- 
sities). These "pertubed" marginal distributions will in general be different 
to the empirical ones obtained from the market. However, the closer these 
"pertubed" marginal distributions are to the empirical ones, the better their 
co-dependence will be modelled by the copula. 

The general formulation of the perturbed copula has been simplified so 
that the whole information of the co-dependence of two random variables 
is condensed in just five parameters with a very intuitive interpretation. 
The skew information of each random variable is introduced through two 
parameters: one controls the volatility level and the other the slope of the 
volatility with respect to strike. Finally, the remaining parameter related to 
the co-dependence is the traditional correlation used for the gaussian copula. 
The current formulation of the perturbed copula only reflects properly the 
skew effect (the slope of the volatility with respect to the strike). In order 
to incorporate the smile effect (slope and convexity) it would be necessary 
to continue the asymptotic expansion up to order two. That's why only two 
parameters (level and slope) are added to each of the random variables in 
addition to the tradicional correlation parameter. 

The contribution of this work is the generalization of the original per- 
turbed copula [1] to cope with lognornal-inspired underlyings, the simplifica- 
tion of the "perturbed copula" general formulation in terms of only five easy- 
to-interpret parameters, their calibration to market data, the application to 
a concrete case (FX quanto options to a third currency), the interpretation of 
the action of the perturbed copula compared to the gaussian copula and the 
comparison of both models with some market standard such the local volatil- 
ity model for a set of scenarios varying market data, moneyness, maturity 
and correlation. 

The organization of this paper starts with the generalization of the origi- 
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nal perturbed copula approach to cope with lognormal-inspired underlyings 
(section [2]) . Section E] reduces the general formulation through some addi- 
tional hypothesis to only depend on five easy-to-interpret parameters. Sec- 
tion |4] shows the calibration procedure of the perturbed copula parameters to 
market data based on a regular Newton-Raphson. Initial parameters based 
on analytical formulas based on approximations calculated through asymp- 
totic expansions are derived in appendix |X1 The interpretation of the action 
of the perturbed copula compared to the gaussian copula is presented in 
section [51 To justify the interpretation, both copulas and a Monte Carlo 
method with local volatility and constant correlation are compared among 
each other through a collection of scenarios pricing FX quanto options to a 
third currency varying skew and correlation. Section E] applies the same tests 
to a real market scenario. Section [7] concludes. 



2 Perturbed copula formulation 

Consider the process {Sl^\ Sl'^\Yt) which follows the dynamics of equation 
(II]), where wi^\ and wP are standard Brownian motions correlated 

with correlations given by equation (|2]), af^ (i=l,2) are the drifts of Sl:^\ m 
is the long term value of Ft, z/ is a parameter which controls the volatility of 
the process Yt and e is a small constant [e « 1) which is the inverse of the 
mean reversion speed (the smaller e, the faster mean reversion). 

dS['^ = ai'^sPdt + Sl'^h {Yt) dWl'^ (1) 
dYt = \{m- Yt) dt + '^dwP 



d{wi^\wl:^^) = pdt d{w^^\wP) = p^Ydt d{wj:^\wj:^^) = p2Ydt (2) 

The processes S'^^^ and Sf^ represent two correlated lognormal-inspired 
underyings. The functions fiiYt) are the volatilities of the underlyings. They 
depend on a common stochastic volatility factor Y^. The correlations p^y 
control the slope of the skew of each underlying and the parameter v controls 
the convexity of the smile. Doing the change of variables X^*'' = hi{S[^^) gives 
the usual normal- inspired processes given by equation ([3]). 
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dXf ) = U''> - Iff (Yt)] dt + h iXt) dwi'^ 

dXp) = [af^ - i/l iYt)) dt + h iXt) dWl'^ (3) 

dYt = \{m-Yt)dt+'-^dwP 

The process Yt is an Ornstein-Uhlenbeck process which is ergodic and 
therefore it has a stationary probabihty distribution (see section 3.2.3 of [2] 
for more information). This means that there exists a probabihty distribu- 
tion so that the expectation of an arbitrary function g of Yt does not depend 
on time (E [(7(1^4)]) = 0). This stationary distribution is reached after 
some amount of time depending on the parameters of the process. For the 
Ornstein-Uhlenbeck process Yt of equation ([3]), this invariant distribution is 
N{m, z/^) (normal with mean m and standard deviation u). See that chang- 
ing the parameter e is equivalent to do a time scale change (the drift term 
is multiplied by 1/e and the diffusion term by the square root of time). In- 
creasing e is equivalent to compress the time scale and thus the result is a 
faster mean reversion. 

The objective is to derive the joint and marginal transition probability 
distributions u^, wf and ti| from an initial point x = {xi,X2) to the end 
point {C,i,C,2), considered constant in this development. These transition 
probability distributions are given by equation (gD, where = (X^^^X^). 
These probability distributions are given by the dynamics of equation ([3]). 

u' = P [x^'^ e da, x^^^ e d6|Xt = ^,Yt = y) 

vl = P (x!^^ E rf^il^f ^ = x^, Yt = y) (4) 

vl = P lx^''> G rf^l^P^ = X2, Yt = y) 

Applying Ito's Lemma to a function which depends on t, x'l^\ x[^^ and 
Yt and setting the drift term to zero (to impose that the evolution of the 
function is a martingale) yields the Kolmogorov backward equation given 
by ([5]), where the operator is given by equation ([6]) in terms of Cq (the 
infinitesimal generator of the Yt process), £1 (the crossed correlation terms 
of processes xl^"* and x[^^ with Yt) and £2 (the infinitesimal generator of the 
joint evolution of xl^"^ and xl^''). These operators are given by equations (I7l) 
to ([9]). The terminal condition for the joint probatility density function u'^ is 
a Dirac delta function 5{^i] xi) with the spike at = Xi. 
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u^{T,xi,X2,y) = 5{^i;xi)6{^2;x2] 



(5) 



= -Co + + £2 (6) 

^o = (m-y)| + .2|, (7) 

£1 = UV^p,yf,(y)^^ + y^p^yf^^y)^^ (8) 

oxioy ox2oy 



£2 = 1 + \n (y) + I/I (y) + Ph (y) h (y) ^ 
+ - i/f (^)) afr + - l/l (^)) ' 



(9) 



dxo 



The solution of equation (E]) is now expanded in powers of y/e according 
to equation (fTOj) . The approximation of the solution will only consider up 
to order 1 {u^ = uq + y/eui). If the solution (fTOj) is substituted in equation 
(|5]) and the terms of the same order are grouped (terms multiplying i, 
1, y/e and so on), equation f|TT]) us obtained. If each term is set to zero 
independently, the system of partial differential equations of equation f lT2|) is 
obtained. This system is easier to solve than equation ([5]). 



u 



' -uo + Veui + eu2 + e^^^ws + ■ ■ ■ (10) 



+ ^ (£oMl + £lMo) + (£oM2 + £lMl + Amo) + .-j^-j^s 
+ (£oM3 + ^lM2 + £2^1) + ■ ■ ■ = ^ 

£o^io = 
LqUx + £iMo = 
£0^*2 + -^iMi + = 
£0^3 + '^lM2 + £2""! = 

As the operator £0 of the first equation of the system (£oMo = 0) is ergodic 
(it is the operator associated with the Ornstein-Uhlenbeck process) and acts 
on y, the solution is constant with respect to y, Uq = Uo(t, xi, X2) (see section 
3.2.4 of [2]). From the second equation of the system f lT2|) (£oMi + £iUo = 0), 
as £1 takes derivatives with respect to y and uq does not depend on y, it 
implies that CiUq = 0. Therefore, £oMi = 0. As again £0 is ergodic and acts 
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on Ui must be constant with respect to y {ui = Ui{t, Xi, X2)). This imphes 
that the combination of the first two terms (^0 + y/^Ui) is independent of the 
volatihty factor y. 

The third equation of the system f|T2|) . reduces to equation f lT3|) as CiUi = 
0, because £1 takes derivatives on y and Ui does not depend on y. This last 
equation is a Poisson equation with respect to the operator £0 in the variable 
y. The general form of this equation is given by ( |T4l) . For this equation to 
have a solution, the operator must be ergodic (with an invariant distribu- 
tion) and the function g{y) must be centered with respect to the invariant 
distribution of the operator. The invariant distribution of the operator Cq 
is a normal distribution A^(m, z/^) and the centering condition reduces to 
(g) = 0, where the notation (■) means the expectation with respect to the 
invariant distribution as given by equation f|T5|) . 

C0U2 + C2U0 = (13) 



Cox + giy) = o (14) 

The centering condition for equation f|T3l) . {C2U0) = 0, gives the zero 
order term of the solution uq. As it does not depend on y, uq can be taken 
out of the expectation. Therefore, uo must satisfy the partial differential 
equation flTBl) . 

{C2)uo{t,Xi,X2) = .^g. 

uo{T,xi,X2) = 6 {^i;x 1)6(^2] X2) 

Taking the expectation (■) over the operator £2 means calculating the 
expectation of each of the coefficients of its partial derivative components as 
they do not take derivatives with respect to y, the variable over which the 
expectation is computed. Equation f|T7|) defines the effective volatilities d"i 
and d"2 and the effective correlation p which appear when this expectation 
is calculated. If the effective parameters are substituted in equation flTB]) . 
equation (fT8|) is obtained. This last equation depends only on the constant 
effective parameters and it is well known. The solution is the transition joint 
density function of two correlated brownian motions, scaled by the effective 
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volatilities and with an instantaneous correlation equal to p. The solution 
is given by equation (IT^ where is given by equation (120]) . It gives the 
transition density of starting at xi and X2 at t and ending at and .^2 at T. 
The volatilities of the two underlying brownian motions are now constant and 
equal to the effective volatilities and the instantaneous correlation is constant 
and equal to the effective correlation. See that the zero order term does not 
depend at all on the stochastic factor Yt. It only provides the contribution of 
the average volatility of each random variable (the effective volatilities) and 
the average correlation. 



0-2 



(/I) 



dug 

at 



— — — 9^ HQ 
' 1 ^ dxidx2 



+ 



2 

E 



Mo (T, xi, X2) = 5 (^1; xi) 5 (^2; X2) 



P 



fa« _ i^2\ duo 



P (/1/2) 

aia-2 



1 -292^0 



(17) 
(18) 



+ j\fds-hj}{T-t) (20) 

As the centering condition (£2^0) = is satisfied, £2^0 is given by equa- 
tion (|2T|) using the definition of the operator ([9]) and taking into account 
that uq does not depend on y (the variable over which the expectation is 
calculated). 



C2U, = C2U0 - {C2U0) = \ iff - (/?)) + I (/I - (/!)) (21) 
+ P (/1/2 - (/1/2)) afj^; - I (/f - (/i )) 1^7 - I (/I - (/I)) ^ 

From equation ( fT3|) the second order order term M2 is given by equation 
(122|), where £0 ^ is the inverse operator (0 = £ ^ (iJ) £0 = H) taken on 
the centered quantity £2 — (£2) given by equation (l2Til . 

M2(t, xi, X2) = -£o ' (£2 - (£2)) uo (22) 

The last equation of the system f|T2|) is again a Poisson equation whose 
centering condition is (£2Mi + £iM2) = 0. Replacing U2 by equation (1221) 
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in the latter centering condition and taking into account that Uq does not 
depend on y yields the partial differential equation (125]) which gives the first 
order term of the solution. The initial condition imposed on this term is zero 
because the initial condition of the solution was already fitted by the zero 
order term. 

(£2) ui{t, Xi,X2) = {CiCq^ {C2 - {C2)))uo ^23) 

Ml(T, Xi,X2) = 

Consider now the functions 4>i{y), 02 (y) and (puiv) as the solutions of 
the Poisson equations fl24l) . See that they are already centered and therefore 
they have a solution. These functions are defined up to an additive constant 
c{t, xi,X2) with respect to y. This constant will be eliminated after applying 
the operator Ci later on in equation fl2B]) as this operator takes derivatives 
with respect to y. 

CoMy) = my) - iff) 

CoHy) = my) - (/I) (24) 
CoMy) = fiiy)f2iy)-{fif2) 

Applying the inverse operator C^^ to equation fl2T|) and taking into ac- 
count the solutions of the Poisson equations (^^, yields equation (^^. If 
the operator £1 is now applied to equation (^^, equation fl2Bl) is obtained, 
where A = CoH^2 - (£2)). 



CoH^2-{C2))uo 



(y) 



iy)t. 



+ 



(y) ^ 



+ 



92 



Mo 



- '2^2 (y) fe 



(25) 



CiAuq = 

+U^p,yh + 0'2a£t| + Wi^^t; - '^'iSf - <^2afe) (26) 

P2y/2 (0;4t; + 0'2|f + 2p0;2afSf - '^lafe - '^'^Sf ) 

Taking expectations of CiA with respect to y to obtain the right hand 
side of equation (l23l) and multiplying by y/e yields equation fl27|) . where 
the constants i?2, -R125 -R215 Q12 and Q21 are given by equation fl28l) 
(see that they all are of the order of y/s). This constants are small and 
will be calibrated from market data (explicit expresions of /«, (f)[ will not be 
necessary) . 
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+ ^lajj^ + R2i^T^ - (Qi2 + Q21) afj^ 

R2 = ^ (/20'2) 

i?12 = (/l0'2) + i^^PP2y (/20'i2) 

i?21 = (/20'l) + iyV2^PPlY {M,2) 

Q2I = ^ {f2<P[) 



(27) 



(28) 



72 

Taking into account the definition of (£2) and tliat equation ( l29l) is satis- 
fied (j + = n), it can be cliecked directly tliat tlie solution of equation ( l23il 
is given by equation (!30|) and therefore, the first order correction is given by 
equation (13T1) . 

ui = -{T -t){CiA)uo (30) 



v^^i = - 7 ^) {^^ (Sf - ) + ^2 (Sf - Sf) (31) 

It is not guaranteed that the asymptotic approximation for the joint den- 
sity function will neither be positive nor integrate 1 in the whole domain. 
To get positive densities the multiplicative expression uq (1 + ta.nh.{y/eui)) is 
used instead of uq + \feu\. The first order Taylor expansions of both expres- 
sions are matched leading to = "^o and u\ = ui/uq. This procedure is 
explained in [1] and it is also shown that the multiplicative approximation 
has the same accuracy than the additive approximation. Therefore, the fi- 
nal solution for the joint probability density function is given by equation 
( 132|) . where uq and y/eui are respectively given by equations (fT9|) and ( 13T|) 
and W is the normalizing constant so that the density integrates 1 in the 
whole domain. In the context of [1], the multiplicative approximation also 
guarantees that the density integrates 1 in the whole domain. Unfortunately, 
the presence of the second partial derivatives (even functions) in ( I3T1) makes 



10 



the density f l32p not integrate 1 and the normahzing constant W has to be 
introduced. Exphcit expressions for the third partial derivatives of Uq can be 
found in appendix A of [1]. 

u'{t, xi, X2; T, ei, 6) = ^uo |l + tanh | (32) 

For the marginal transition density functions vf and v^, the same argu- 
ment used to obtain the solution for the joint distribution can be applied. 
The solution is given by equation (!33|) where the zero order approximation is 
a regular normal density function given by equation where Xj is again 
given by fl20l) and Wi is the normalizing constant so that the marginal density 
integrates 1 in the whole domain. 



W, 



Pi 



1 + tanh -(T-t) 



Ri (d^pi d'^pi 



Pi 



(33) 



Piity ^ii T, 



(JiJ2Ti{T-t) 



: exp 



{ii - ^ 



~ \2 



(34) 



The perturbed copula function will be given by the ratio of the joint 
density and the product of the two marginal densities as shown by equation 
where zi and Z2 are related to ^1 and ^2 by equation The variables 
Zi (they are in the interval [0, 1]) are the marginal cummulative probabilities 
of ^j. The copula function is expressed in terms of Zi and for each of those 
values, the values of should be calculated through the quantile function 
(the inverse of the cummulative density function). It is proved in that 
fcop is indeed a copula function. 



fcop {zi,Z2) 



U%t,Xi,X2]T, ^1,^2) 

vl{t,XuT,^i)v^{t,X2,T, 5) 



p [x!p < e 



y 



(35) 



(36) 



See that for the purpose of calculating the copula density, it can be as- 
sumed that Xi = 0, because the quantile function provides a value = ^i+Xi, 
where C,i is the increment of C,i with respect to the mean x,. When this value 
is replaced in the terms (^j — Xj) of equations (fT9|) and the result is 

{C,i — Xi) = ^i + Xi — Xi = ^i. Therefore, if we had started with instead of 
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the result would have been the same (for numerical reasons it is better to 
assume Xj = for the purpose of calculating the copula function). 



3 Reduction to five intuitive parameters 

This section incorporates some hypothesis in order to express the perturbed 
copula in terms of just five easy-to-calibrate parameters with a clear in- 
terpretation. Equation f l37|) shows the main hypothesis for this parameter 
reduction: the volatility dynamics of both random variables have the same 
dependence with respect to the common volatility factor Yf, only the level of 
volatility might be different. This common dependence with respect to Yf is 
expressed through the function g{y). 



fiiy) = (^i9{y) f2iy) = (y2g{y) (37) 

If the functions /j from equation (|37|) are replaced in equation (|T7|) and 
the new hypothesis {g"^) = 1 is added, the expression for ai and p turns very 
simple into equations ( l38l) and ( l39l) . 

^i = \/(/f)=^i\/(^ = ^i ^2 = V^=^2V^ = cr2 (38) 



P = — — = P 1— 1— = P (39) 



Using the linearity property of the operators, the solutions of the three 
Poisson equations i^^, 0i, 02 and 0i2 can be expressed in terms of a single 
solution (f) according to equation ( HOl) which satisfies the Poisson equation 

m- 

My) = ^10(2/) My) = ^IM) M{y) = (^1(^2^) (40) 



CoM) = 9\y) - (/> (41) 

Then the parameters i?2, -^12, -R215 Q\2 and Q2\ from equation f l2S]) 
turn into equation P2|) . where i?2, c"i and will be calibrated to market 
(see section H]) and i?i2, -R21, Q12, Q2\ will be calculated according to equation 

(HID. 
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R2 = {g4>') 

R12 = (^^criai + uV2^pp2Ycrla,) {g<P') = (^)^ + 2 R,p 

R21 = (^^2^? + uV2^pp,Yala,) {g<P') = l^yR, + 2 (f^) R^p ^^^^ 

Q21 = '-^cr,a! {g<P') = R, 

Section H] along with appendix |X] will show that the marginal densities 
given by equation fl33|) produce an implied volatility curve that is a straight 
line with respect to a given definition of moneyness. This line is defined by 
the slope and the intercept (two parameters). Therefore, the parameters (Tj 
and Ri {i = 1, 2) govern respectively the ATM (at-the-money) level and the 
skew slope for each of the marginal distributions. The parameter p mea- 
sures the correlation of the co-dependence (it has the same meaning for the 
gaussian copula). The parameters ai = ai govern the volatility level because 
they are the volatilities for both the joint and marginal density functions of 
equations (IT^ and provided by the zero order term of the asymptotic 
approximation. The parameters Ri govern the slope of the skew because for a 
fixed level of (jj, they only depend on the correlation piy and some common 
parameters {v, {Q'P') and \/s). The correlation p^y between the brownian 
motions driving the two random variables Wt^'^ and the common factor Y-t 
controls the slope of the skew. Finally, the parameter p controls the corre- 
lation between both random variables because it is the correlation for the 
zero order term of the joint density function according to equation fll9p . The 
hypothesis considered assume that the internal dynamics for both random 
variables are the same. However, the side and amount of skew is independent 
for each variable (it is controlled independently by each correlation Piv)- 



4 Calibration 

The calibration procedure involves finding the parameters cxj (implied volatil- 
ity level) and Ri (implied volatility slope) for both underlyings (2 degrees of 
freedom for each underlying). These values along with the correlation p 
(which is not calibrated but input by the trader) will allow calculating the 
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joint and marginal densities and therefore the copula function. The method 
proposed is an exact calibration to two vanilla options for each underlying 
using a simple Newton-Raphson algorithm. However, for the algorithm to 
converge, initial values close enough to the solution are needed. These initial 
values are calculated using asymptotic expansions for vanilla option prices 
similar to those described in section H] for the pertubed copula. These initial 
values are given by a calibration procedure described in [2], but this cali- 
bration is carried out for processes formulated in the real market measure 
(rather than the risk free measure), assuming a risk premium different from 
zero. Therefore, in order to better relate the calibration procedure of [2] with 
the estimation of the initial parameters, this section formulates the process 
in the real measure. 

Equation fH3l) presents the evolution of the underlying processes S^^ (for 
i = 1,2), where /if"* is the unknown market drift of the underlying, gf"* is 
the continuously compounded dividend yield and the rest of the parameters 
are the same as those defined in equation ([1]). The asterisk in the brownian 
motions expresses that process is referred the real market measure rather 
than the risk neutral measure. The noise driving the stochastic volatility 
dW*^"^^ has already been decomposed into a linear combination of indepen- 
dent brownian motions dW^ and dZ^ to fix the appropriate correlation piY 
between the underlying and the stochastic volatility brownian motions. 

dSj^^ = (/xf ) - g«)5«rft + f^iY^sj^^dW: 

dYt = l{m-Yt)dt + ^ {^PiYdW: + ^1 - plydZ;^ ^^^^ 

dWt = dW; + ^',^1' dZt = dZ; + -^dt (44) 

Equation (jH]) shows the change of measure to turn the real market mea- 
sure into the risk neutral measure. See that ^j(y^)''' is the market risk pre- 
mium and 7j is the risk premium for volatility risk. When dW* and dZ^ are 
replaced in equation fH5]) . equation fH5]) is obtained. It expresses the evo- 
lution of the underlying in the risk neutral measure. The function A(yj) in 
equation ( H6|l is a combined market premium and volatility risk. See that the 
unknown drift pf^ of Sf'' turns into the known risk free rate. Therefore, the 
process of the discounted underlying (S'^*-' exp(— /q r^fis)) will be a martin- 
gale (zero drift) which is the condition required for the absence of arbitrage. 
When only pricing is considered, the risk premium is set to zero (A(yf) = 0). 
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dS^ = in - qi^)Sfht + m)Sf^dW, 

dYt ={^{m- Yt) dt - (Yt)) + ^ [p^ydWt + y'l - pfydZ^j ^^^^ 

A {y) = pt ^f^Y,)' + ^^^^ 

The calibration of the process of equation fH5]) is described in |2]. How- 
ever, it was done considering that the process was not written in the risk 
free measure but the real market measure. Indeed, [2] estimates the effec- 
tive volatilities cij of equation ( !T7|) from historical returns whereas section [3] 
sets them equal to cij according to equation (155]) under the assumption that 
{g"^) = 1. This can be done because in the context of this paper, the pro- 
cess will be only used in the risk neutral measure (the market risk premium 
will be set equal to zero) and therefore, no historical data will be needed to 
estimate the level and slope of the implied volatility. 



Ri = -aaf = ai^h- — 48 

a a I 

Equation (H7|) shows the approximation up to first order of the implied 
volatility al^^ of a vanilla option (see appendix |A] to see where this expres- 
sion comes from) when the underlying follows the process of equation (1451) 
assuming that e is small (fast mean reversion). Fj^ = Si^ exp(//' (x§ ds) is the 
forward value of the underlying i at time T. For pricing purposes, the risk 
neutral meaure will be used and a^*-* = (r^ — g^) where is the domestic risk 
free rate and is the dividend yield (or the foreign risk free rate if an FX pair 
is considered). See that up to first order, the implied volatility behaves as a 
straight line where a is the slope and h is the intercept, when the independent 
variable is equal to the log-moneyness-to-maturity ratio (In (;^') Y^t>- "^^^ 
parameters a and h are estimated through a linear regression of the volatility 
with respect to the log-moneyness-to-maturity ratio. Equation fHHl) gives the 
parameters of the model Ri and Oi in terms of a and h. The second order 
equation for Oi can be simplified assuming that Oi ~ 6, yielding the right 
expression of equation fHS]) . 
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poo 

Ci{K, T) = P{t, T) / (exp (A + e.) - Ktvl (t, 0, T, ^i) d^i (49) 



Ao = In + ^ (r» - gW)^^. - -erf (T - t)^ (50) 

Equation (l49ll shows how to calculate the price of a call option out of 
the perturbed marginal transition probability density v| given by equations 
(IMl) and ( l33l) . where P{t,T) is the discount factor from the expiration date 
T to present time t and /3j is a parameter to fix the mean of the distribution. 
See that the transition probability is defined from to (it has no mean), 
because the mean of the distribution is taken into account by Pi. An initial 
value for this parameter /3jo is given by equation f lSU]) which is the mean of the 
zero order term of the solution given by equation f l20p . where af ^ = r^*'' — g^*'* 
and Xi = ln{St^^). 

The parameter /3jo is only an asymptotic approximation of the actual pa- 
rameter Pi which satisfies equation fl5T]) (the expected value of the underlying 
must be equal to its forward value). The parameter Pi is easily obtained cal- 
culating the forward from the initial Pio and integrating numerically 
equation ( 15T1) . This value will not in general match the actual forward FiT. 
However, if the next beta is corrected by /Sf*^^* = /Jf^" — ln{F^^^) + In^Fix) 
the new value of the forward F^'^^ obtained out of P'^^^^ will be a lot closer 
to FiT- After a few iterations. Pi converges. 

/oo 
exp(A + e.)<(t,0,T,6)c?6 (51) 
-oo 

Once the forward of the distribution FiT is matched, the price of two 
vanilla options can be matched using a simple Newton-Raphson method with 
step shortening starting from the initial Ri and (jj given by equation (HHll . 

5 Interpretation of perturbed copula 

This section interprets the calibration and the effect of the perturbed copula 
compared to the gaussian copula. In terms of pricing, the interpretation is 
carried out applying the perturbed gaussian copula approach to the valua- 
tion of FX quanto options to a third currency different to the currencies of 
the underlying pair. In particular, the underlying FX pair considered is the 
XAU/USD which is quoted in USD per ounce of gold (XAU). The option is 
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quantoed to EUR and therefore the second underlying pair involved is the 
EUR/USD. Equation ( 15^ shows the payoff function of this option, where St 
is the price of the XAU/USD (in USD per ounce of gold), Xt is the price of 
the EUR/USD (in USD per unit of EUR), DF^^^ is the discount factor of 
the USD curve and K is the strike price. If the USD money market account 
is chosen as numeraire, both 5*^ and are denominated in the numeraire 
currency (USD) and their drifts are simply calculated as the difference be- 
tween their domestic (USD) and foreign interest rates at maturity. However, 
as the option is quanto, the {St — K)~^ will be paid in EUR and therefore, 
as it should be discounted with the numeraire in USD, the quantity must 
be converted to USD first, multiplying by the EUR/USD exchange rate Xt- 
The spot price of the XAU/USD is Sq = 981.3 and the spot price of the 
EUR/USD is Xo = 1.422. 

p = E [{St -KY XtDF^^'^] (52) 

The expectation given by equation (152|) is computed through the double 
integral given by equation (153|) . where are the logarithms of the underlyings 
at expiration (^i = In(S'T) and ^2 = ^'^{Xt)) and /J"*"-* is the joint probability 
density function of both underlyings given by equation (see chapter 
2 of to see where this equation comes from), where fcop is the copula 
function defined in equation fl35|) . zi are given by equation fl36l) and //""'^^ 
are the empirical marginal distributions. The latter are obtained through 
equation flS^ (see |1] to find out where this equation comes from), where 
Pi{K, T, cr*™'P') is the Black Scholes price of a put option of underlying i with 
strike K, maturity T and al"^^\K,T) is the interpolated implied volatility 
for strike level K of underlying i at maturity T. 



p = DF^'^ J (e«^ - KYe^-p^'^' (6, 6) d^id^2 (53) 

R2 

/^™°* (6, 6) = fcop {zu Z2) fr' (6) ir' (6) (54) 

The perturbed copula has been tested in a set of 25 scenarios varying 
the correlation and the skew of both underlyings. These scenarios have zero 
interest rates and have been created with volatility surfaces given out of a 
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Heston model with constant parameters varying the correlation between the 
brownian motions of the underlying and the variance process (this parameter 
controls the skew). The time horizon considered is one year. A skew which 
favors out-of-the money (OTM) puts (lower strikes have higher volatility) is 
referred to in this paper as left skew. When OTM calls are favored (higher 
strikes have greater volatility), the resulting skew is called right skew. There 
are five groups of scenarios with different skews: "LR" (left-right), "RL" 
(right-left), "SS" (smile-smile), "RR" (right-right) and "LL" (left-left). Each 
group has five differerent correlations between both underlyings (0.6, 0.3, 0, 
-0.3, -0.6). The perturbed copula has been calibrated to the 25 delta out-of- 
the- money call and put options (the most liquid products). 
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Figure 1: Implied volatility of left and right skew Heston scenarios compared 
to their calibrations (left plots), corresponding perturbed pdf (middle plots) 
and empirical Heston pdf (right plots). Calibrated parameters: /3i = 6.8783, 
(Ti = 0.0892, Ri = 1.31 ■ 10"^ (up) and /^a = 0.3550, (Xa = 0.0877, Ri = 
-1.29 ■ 10"^ (down). 

Figure [H presents the left and right skew scenarios. The upper plots show 
the left skew scenario (higher volatility for lower strikes) and the lower plots 
show the right skew scenario (higher volatility for higher strikes). The cali- 
brated parameters for the left skew scenarios are f3i = 6.8783, Ui = 0.0892, 
Ri = 1.31 ■ lO"'' and for the right skew scenario are (32 = 0.3550, (Xa = 0.0877, 
i?2 = — 1.29 ■ 10"*^. See that the f3i that appears in equation (149|) is a lot 
bigger than /^a because this parameter fits the mean and the spot levels are 
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quite different (981.3 versus 1.422). On the other hand, the volatihty level is 
not very high (around 9%), the skew levels are rather mild and have opposite 
signs (in agreement with smirks opposite to each other). The left plots of 
figure [1] show the implied volatility of the original Heston-generated (labeled 
"Original" ) and the perturbed-copula (labelled "Perturbed" ) calibrated sur- 
faces versus the log-moneyness-to-maturity ratio {hi[K/FT)/T, where Ft is 
the forward of the underlying at maturity T) on the horizontal axis. The mid- 
dle plots show the perturbed marginal probability density functions obtained 
after cahbration from equation ( 133|) (labeled "Pertubed" ) and the zero order 
term from equation (labeled "Lognormal" and centered around /3j) also 
displayed for comparison purposes (these functions can be replicated with 
the parameters mentioned at the beginning of the paragraph). The right 
plots show the empirical marginal density functions obtained from equation 
fl55|) out of the Heston-generated volatility surface. See that marginal den- 
sities are expressed in terms of the logarighm of the underlying rather than 
the underlying level itself. It can be seen that the fit of the skew provided 
by the perturbed marginal density function is very reasonable. It is clear to 
interpret that the left skew scenario increases the weight of the left queue 
(lower underlying values) and displaces the mode of the distribution to the 
right (higher underlying values). The opposite happens for the right skew 
scenario (increase of right queue density and move of the mode of the dis- 
tribution to the left). See that the perturbed density functions incorporate 
a slight bump in the queues to create the skew. For extreme queue values, 
the density function is not as fat as the that of the empirical distribution 
and that is why the skew provided by the perturbed copula flattens for very 
extreme values of the queue (this will be clearer in section [6]). 

From the five groups of scenarios and 5 correlations, only one of them (the 
"LL" scenario) for the extreme correlation values (0.6 and -0.6) is analyzed. 
The rest of combinations are symmetrical with respect to this reference sce- 
nario. The left plots of figure [2] show the perturbed joint probability density 
function (numerator of the copula function) of both underlyings given by 
equation (152]) . The middle plots of figure [2] present the copula function given 
by equation (1351) for values greater than 1 (to allow for clarity). Finally, the 
right plots of figure [2] display the ratio between the pertubed and gaussian 
copula functions also for values greater than 1. These plots might be repli- 
cated using the calibrated parameters given in figure [1] and calculating Rij 
and Qij according to equation |l2l The upper plots correspond to the positive 
(0.6) correlation scenario and the lower plots to the negative (-0.6) correla- 
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Figure 2: Copula numerator (left), copula function greater than 1 (middle) 
and ratio greater than 1 between perturbed and gaussian copulas (right) for a 
LL skew Heston scenario with positively (upper plots) and negatively (lower 
plots) correlated underlyings. 




Figure 3: Joint pdf (copula function times empirical marginal densities) ob- 
tained with gaussian copula (left), perturbed copula (middle) and ratio be- 
tween them greater than one (right) for a LL skew Heston scenario with 
positively (upper plots) and negatively (lower plots) correlated underlyings. 
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tion scenario. The upper left plot of figure [2] shows a positively sloped (from 
bottom left to up right) oval form which corresponds to a bi-normal distribu- 
tion with positive correlation. However, the mode of the distribution is not 
situated in the center but it is displaced upwards to the right. A "LL" sce- 
nario increases the left queues (lower values) and shifts the mode to the right 
(higher values) of both underlyings. It can be seen that the displacement 
produced by the skew goes to the upper right corner in the same direction 
as the correlation (the queues are displaced in the opposite direction towards 
the lower left corner). For the negative correlation scenario (bottom left plot) 
the displacement of the skew goes perpendicular to the direction of the cor- 
relation. The middle plots of figure [2] show the copula function. When the 
copula function is equal to one, the joint probability density function is the 
product of the two marginals indicating independence or no co-dependence. 
When the copula function is greater than one, the density is increased indi- 
cating that there is more co-dependence (the opposite happens with a copula 
function less than one). 

The middle upper plot of figure [2] shows how the copula function increases 
the probability in the main diagonal as the correlation is positive (it is not 
clear in the plot but the lower left corner increases the density more than the 
oppossite corner for the effect of the skew) . The middle lower plot of figure [2] 
shows that the copula is greater than 1 in the anti-diagonal as the underlyings 
are negatively correlated. However, the anti-diagonal is displaced or biased 
through the the lower left corner because the effect of the "LL" skew scenario 
is moving the queues in that direction. 

The right plots of figure [2] show the ratio between the perturbed and 
gaussian copulas and allows the comparison between them. The upper right 
plot corresponds to the positive correlation and shows that the effect of the 
perturbed copula is to increase the density of the left queues of both under- 
lyings. For the horizontal variable, the left queue density is increased for 
the points with less co-dependence (the upper left points increase the most 
and progressively increase less until the lower left corner is reached) . For the 
vertical variable, the left queue density (the lower side of the plot) is again 
increased more for the points with less co-dependence (lower right side of the 
plot) and progressively decreases as the points with more co-dependence are 
reached (lower left corner). The lower right plot of figure [2] corresponds to the 
negative correlation scenario. Again, the left queue densities are increased as 
it can be seen with the significant increase of density in the lower left corner 
of the plot, precisely where the co-dependence is the smallest. For higher 
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values of the co-dependence (the anti-diagonal), the densities are displaced 
towards the lower left corner. 

Figure [3] compares the joint probability functions given by equation fj54|) 
using the gaussian copula (left plots), the perturbed copula (middle plots) 
and the ratio of both (right plots) for values greater than zero (plots are 
clearer this way). The upper plots correspond to a scenario with positive 
correlation whereas the lower plots correspond to negative correlation. The 
upper left and middle plots of figure [3] shows the positive correlation oval 
shape through the diagonal. The distribution obtained with the pertubed 
copula shows slightly displaced or deformed towards the upper right corner 
(as corresponds to the movement of the mode of the joint distribution as seen 
in the upper left plot of figure [2]) . The upper right plot of figure [3] shows how 
the left queues in the region with less co-dependence have more density for 
the perturbed copula than the gaussian copula (this is the same as happened 
with the upper right plot of figure [2l The lower left and middle plots of 
figure |3] show the anti-diagonal oval shape which corresponds to negatively 
correlated underlyings. Again, the distribution obtained with the perturbed 
copula appears as a deformation of the gaussian copula towards the upper 
right corner. On the other hand the lower right plot of figure [3] shows that the 
density is increased in two regions aligned with the anti-diagonal direction. 
Looking at the lower right plot of figure [2] (the ratio of the perturbed and 
gaussian copula functions), the mentioned regions correspond to the density 
increase of the perturbed copula in the anti-diagonal belt. 

The corresponding plots given by figures [2] and [3] for the "LL" scenario 
which would result for the rest of scenarios can be obtained through sym- 
metries. For a "RR" scenario, the right queues of both underlyings would 
be fatter going to the upper right corner and the mode of the joint distri- 
bution would be deformed in the opposite direction (towards the lower left 
corner) and the plots of figure [2] would be symmetrical with respect to the 
anti-diagonal line. The plots of a "RL" scenario with positive correlation 
would be symmetrical to the lower plots of figure |2] with respect to a vertical 
axis situated in the middle of the plots (the mode of the distribution would 
be deformed towards the upper left corner). The same symmetries apply to a 
negative correlation scenario but the plots would be symmetric to the upper 
plots of figure |2] with respect to the same axis. The plots of a "LR" scenario 
with positive correlation would be symmetrical to the lower plots of figure 
|2]with respect to a horizontal axis situated in the middle of each plot (the 
joint distribution would be deformed towards the lower right corner and the 
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queues to the upper left corner) . The same symmetries apply to the negative 
correlation scenario but the plots would be symmetric to the upper plots of 
figure [21 

If the asymptotic approximation of the perturbed copula were carried 
out up to second order, the smile effect could be captured. In this situation, 
the perturbed copula would increase the density around the mode of the 
distribution as well as both queues and will reduce the density in between. 
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Table 1: Comparison of Gaussian and perturbed copula methods. 

Table [1] compares the prices given by the gaussian and perturbed copulas 
of equation (!52|) for the given set of 25 scenarios varying skew configurations 
and correlation. The strike price is K = 834.105 and it is in-the-money. The 
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columns present the prices of the gaussian (" Gcop" ) and perturbed (" Pcop" ) 
copulas and the difference of the perturbed minus the gaussian copulas (" P- 
G") in per unit of notional (e.g. 0.0001 is 0.01% of notional or a basis 
point). The "SS" scenarios (11 to 15) have been omitted because they give 
no correction (the pertubed copula can only capture the skew effect but not 
the smile). The interpretation in terms of pricing is in general not as clear. 
The perturbed copula provides positive corrections for the "LL" scenarios (21 
to 25). This is not surprising as the joint distribution of figure |3] is deformed 
towards the upper right corner (the direction in which the payment increases) 
and the queues don't pay out as it is a call option. For the positive correlation 
case (upper right plot of figure [3]), the queues increase the density in areas 
with higher payoff (the left queue of St does not pay out as a call option is 
being priced). By looking at the negative correlation case (lower right plot of 
figure [3]) it is not as clear from a qualitative point of view that the correction 
should be positive. However, table [T]shows that indeed it is positive and even 
greater than the positive correlation case. The corrections for scenario "RR" 
(the symmetrical of "LL" ) are of opposite sign and about the same magnitude 
as the correponding "LL" scenarios. The "LR" scenarios involve deforming 
the joint distribution towards the right lower corner of the distribution (the 
fatter queues move in the opposite direction towards the upper left corner). 
The left tail of St does not pay out (it is a call option) and therefore the 
mode of the distribution moving to the right increases the payoff and the 
left tail getting fatter does not reduce the price as the payoff is zero. On 
the other hand, the payoff for Xt gets lower on one side as the density is 
greater for lower values but also gets higher on the other as the density is 
increased for very high values (the right tail gets fatter). These effects might 
compensate with each other or even the effect of right tail (higher values 
of Xt) be more significant. This means that it is reasonable to think that 
the perturbed copula gives a positive correction for a "LR" scenario (see 
scenarios 1 to 5 in table [1]). Following a similar reasoning, the symmetric 
scenario "RL" should give opposite corrections (see negative corrections of 
scenarios 6 to 10 in tabled]). 

This set of 25 scenarios was extended to cover three moneyness levels 
with strikes K = 834.10,981.30,1128.50 (in-the-money, at-the-money and 
out-of-the-money) and two maturities (1 and 2 years). This new set of sce- 
narios was compared with a standard local volatility model (see for example 
[3] and |1]) with constant correlation set equal to the correlation used for 
the copulas. This means that the noises of the two underlyings were corre- 
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lated using constant correlation and the distribution of the underlyings was 
obtained simulating each of them with their corresponding local volatility 
until expiration. The comparison showed that the prices obtained with the 
gaussian copula were not any further than 10 basis points out of the whole 
set of scenarios from those of the Monte Carlo method with local volatil- 
ity. For this set of scenarios, the maximum differences between the gaussian 
and perturbed copulas were obtained for the 2 year maturity and the in-the- 
money cases and it was 92 basis points. The conclusion of this study is that 
a standard local volatility model is rather equivalent to the gaussian copula 
and the corrections provided by the perturbed copula can reach up to almost 
100 basis points (1% of the notional). 

6 Case study: FX quanto options 

This section compares the perturbed copula, the gaussian copula and the 
Monte Carlo method with local volatility and constant correlation for a set 
of scenarios build out of a real market scenario for the same FX call option 
on XAU/USD quantoed to EUR considered in section [51 This scenario cor- 
responds to a "LS" scenario with the XAU/USD highly left skewed and the 
EUR/USD very mildly right skewed (almost a smile). Five correlations (0.6, 
0.3, 0, -0.3, -0.6), two maturities (1 and 2 years) and 5 moneyness levels 
(0.7, 0.85, ATM, 1.15 and 1.2) with strikes K = (656.46, 797.13 , 937.79, 
1078.47, 1125.36) were considered. The spot price of the XAU/USD and the 
EUR/USD are So = 937.79 and Xq = 1.4029. 

Figure m shows the 1 year implied volatilities (left plots), the calibrated 
perturbed marginal density (middle plots) and the empirical market marginal 
density (right plots). The upper plots correspond to the XAU/USD pair 
and the lower plots to the EUR/USD. The calibration parameters for the 
XAU/USD are /3i = 6.7723, ai = 0.2989 and Ri = 73.43 ■ 10"^. For the 
EUR/USD they are (32 = 0.3287, ^2 = 0.1492 and Ri = -1.08 ■ 10"^. See 
that the volatility levels are now more realistic than those of section |5] (instead 
of around 9% they are now around 30% and 15%). The skew level of the 
XAU/USD is very high (around 50 times the skew levels of section [S]). The 
left plots of figure m show the calibrated (labelled "PerCop") and the original 
market (labelled "Original") implied volatilities and the middle and right 
plots show the pertubed and empirical density functions (this figure has the 
same structure than figure S]). It can be seen that the calibration of the 
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Figure 4: 1 year implied volatility of left skew and smile market scenar- 
ios compared to their calibrations (left plots), corresponding perturbed pdf 
(middle plots) and empirical market pdf (right plots). Calibrated parame- 
ters: /3i = 6.7723, (Ji = 0.2989 and Ri = 73.43 ■ 10"^ (XAU/USD up) and 
/32 = 0.3287, as = 0.1492 and Ri = -1.08 ■ 10"^ (EUR/USD down). 

perturbed marginal density (see upper left plot of figure H]) of the XAU/USD 
to market is reasonable (very left skewed). If the perturbed and empirical 
densities are compared (upper middle and right plots), it can be seen that 
the pertubed copula skew is rather extreme and shows almost a bimodal 
distribution with a remarcable bump in the left queue (as already mentioned, 
the skew of the XAU/USD is around 50 times the skew considered in section 
|5]). This bimodal distribution fits the skew up to a level of moneyness of 
around -0.4. For lower moneyness levels, the skew flattens as the queue is 
not as fat as the empirical distribution for very low values of the underlying. 
Looking at the lower plots of figure HI the calibration of the EUR/USD is 
rather poor but still mildly right skewed. This is the consequence of the fact 
that the perturbed marginal density can only capture skew but not smile. 
The empirical market density (lower right plot) has fatter queues and higher 
mode point than the perturbed density. However, the perturbed density has 
higher density in between the queues and the mode. 

Figure [5] has the same structure than figure [2] and presents the 1 year 
perturbed joint density functions given by equation (15^ (left plots), the per- 
turbed copula functions given by equation (!35|) (middle plots) and the ratio 
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Figure 5: 1 year copula numerator (left), copula function greater than 1 
(middle) and ratio between perturbed and gaussian copulas (right) for val- 
ues greater than one considering a LS skew market scenario with positively 
(upper plots) and negatively (lower plots) correlated underlyings. 



Mkt1ylTM1:GaussCop joint pdl MktlylTMl: PerCop joint pdf MktlylTMl : Ratio joint pd! PC/GC > 1 




400 600 800 1000 1200 1400 1600 400 600 800 1000 1200 1400 1600 400 600 300 1000 1200 1400 1600 

1 st underlying spot 1 st underlying spot 1 st underlying spot 



Figure 6: Joint pdf (copula function times empirical marginal densities) 
obtained with gaussian copula (left), perturbed copula (middle) and ratio 
between them (right) for values greater than one considering a LS market 
scenario with positively (upper plots) and negatively (lower plots) correlated 
underlyings. 
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Table 2: Perturbed minus gaussian copula difference in basis points varying 
moneyness, correlation and maturity for market scenarios. 
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Table 3: Local volatility with ATM correlation minus gaussian copula differ- 
ence in basis points varying moneyness, correlation and maturity for market 
scenarios. 
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of the perturbed over the gaussian copula for values greater than one (right 
plots). Again, these plots might be replicated using the calibrated parame- 
ters given in figure [Hand calculating Rij and Qij according to equation fH2|) . 
The upper and lower plots correspond respectively to positive and nega- 
tive correlation between the XAU/USD (horizontal axis) and the EUR/USD 
(vertical axis). The upper and lower plots are approximately symmetrical to 
each other with respect to a horizonal axis in the middle of the plots. This 
is because the skew with respect to the EUR/USD (the vertical axis) is very 
mild (almost a smile). The lower plots of figure [5] are very similar to the 
lower plots of figure [2l Now the joint distribution (lower left plot of figure [5]) 
gets deformed only to the right (and not upwards to the right). In addition, 
only the left tail of the horizontal axis is fatter (the lower middle plot has 
higher copula values for the upper left corner and not as big for the opposite 
corner). Again, the lower right plot of figure |5] shows how the left tail den- 
sity increases more for the points with less co-dependence (lower left corner 
has the biggest increase in density with respect to the gaussian copula) and 
progressively increases less as the co-dependence rises. 

Figure [H] presents the final joint density function from equation f l54l) for 1 
year maturity obtained with the gaussian (left plots) and perturbed (middle 
plots) copulas and the ratio of the perturbed over the gaussian density for 
values greater than one (right plots). The upper plots correspond to positive 
correlation and the lower to negative correlation between the XAU/USD and 
the EUR/USD. The effect of the left skew can be seen comparing the left 
and middle plots of figure (6) the distribution gets deformed towards the right 
for the perturbed copula case compared to the gaussian copula (the mode of 
the distribution is moved towards the right and the tail gets fatter towards 
the left). The left tails of the XAU/USD rise the density starting from 
the points with less co-dependence and progressively increasing less when 
approching the points with greater co-dependence (these points also increase 
their density but as they already have big density, the rise by the effect of the 
skew is less than for the points with almost no co-dependence). The upper 
right plot of figure E] shows a big area of increase of density in the upper 
middle side of the plot. This is the effect of the bump which appeared in the 
marginal density function (see the left tail of the pertubed marginal density 
of figure H]). For the negative correlation scenario (lower plots of figure [6]), a 
symmetric behaviour shows up with respect to a horizontal axis through the 
middle of the plots. 

Table [2] shows the difference in basis points between the perturbed and 
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gaussian copulas for the "LS" scenario varying moneyness, maturity and 
correlation. In order to properly compare perturbed and Gaussian copulas, 
the correlations which appear in the columns of table [2] refer to those input 
in the Gaussian copula. The correlations used for the pertubed copula are 
implied so that the quanto forward, E (StXt), is equal for both the Gaussian 
and perturbed copulas. The upper group of rows show different moneyness 
levels for 1 year maturity and the lower group of rows for 2 year maturity. 
Differences increase with maturity and they can get beyond 40 basis points 
(0.4% of the notional). As the overall effect of the skew of the XAU/USD is to 
deform the joint distribution towards higher values of the horizontal axis (in 
this direction the payoff is higher), the corrections given by the perturbed 
copula are positive for almost every case. The fact that the left tail gets 
fatter does not have a big impact as the payoff for the left tail is zero (a call 
option is considered). It is not easy to interpret the premium corrections in 
terms of the correlation (whether positive correlation should correct more or 
less than negative correlation). In fact, the same qualitative plot given by 
figure E] appears for scenarios which invert the size of corrections between 
positive and negative correlations (all corrections were however positive for 
these cases). 

Table |3] is similar to table [2] but shows the difference of the Monte Carlo 
method with local volatility and constant correlation minus the gaussian 
copula. See that the biggest difference does not go beyond 17 basis points. 
This means that the Monte Carlo method is rather equivalent to the gaussian 
copula (as it was concluded in section [5] varying moneyness, maturity and all 
combinations of skew) . 

7 Conclusions 

The perturbed copula approach of [T] has been successfully applied for val- 
uation of derivatives which depend on two underlyings for which prior in- 
formation about their co-dependence is unknown. The application of this 
perturbed copula allows to introduce the skew information of the underly- 
ings in the co-dependence. The perturbed copula formulation has been gen- 
eralized to deal with widely used lognormal-inspired underlyings where the 
skew is introduced through the dependence of the volatility with respect to a 
common stochastic volatility factor. The analytic formula of the perturbed 
copula is obtained through an asymptotic expansion under the assumption 
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that this common factor is fast mean reverting. 

The original formulation of the peturbed copula has been simplified through 
additional hypothesis which allow to reduce the copula to five intuitive and 
easy to interpret parameters: two control the volatility levels, two the skew 
levels and one the correlation between both underlyings. 

An exact fit calibration procedure of this five parameters is proposed 
giving initial values for them which are close enough to the solution so that 
a regular Newton-Raphson search algorithm finds the solution quickly. 

The effect of the perturbed copula is interpreted in comparison with the 
gaussian copula in terms of the direction of the skew and the correlation 
of each underlying. Intuitive criteria are provided to qualitatively predict 
the effects of the perturbed copula in the pricing of a derivative. These 
qualitative ideas are applied to the particular case of FX options quantoed 
to a currency different from the currencies of the underlying pair. 

A real market case study is also analyzed for FX quanto options to a third 
currency. It has been seen that the price impact of considering the skew in 
the co-dependence is not negligible (it can go beyond 40 basis points or 0.4% 
of the notional amount in some of the cases analyzed). This rises a warning 
concern about the risk of quanto options in the presence of skew. 

The cases analyzed have also been compared with a Monte Carlo local 
volatility model with constant correlation equal to the correlation used in 
the copulas. It has been concluded that this method is almost equivalent to 
the gaussian copula. This means that the regular widely used local volatility 
models do not incorporate the skew effect in the co-dependence and some 
potential model risk might be unidentified. 
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A Initial parameters for calibration 

This appendix derives the initial guess of the calibration parameters. This 
derivation has already been carried out in [2| and therefore it will be omitted 
here. However, there are a few subtle differences in the derivation process 
which are addressed in this section. 

Consider the system S'^*^ and Yt of stochastic equations psl) . The same 
procedure used for the derivation of the perturbed copula is followed here to 
obtain the asymptotic expansion of the vanilla option price. The operators 
£i and £2 given by equations ([H]) and are replaced by £1 and £2 given by 
equations (156!) and (!57|) . The operator £0, which corresponds to the Ornstein- 
Uhlenbeck process, does not change. See that the operator £2 is the operator 
Cbs resulting out of the Black Scholes process with volatility equal to f{y)- 
If the underlying process is denoted by Sj:^^ and the underlying variable by 
Si, Si multiplies the partial derivative g^^^ and Sf multiplies ^ because 
the lognormal-inspired process is considered instead of the normal process of 
equation used for the derivation of the copula. The term (r — q)Si-^ 

corresponds to the drift of the process Sl^\ The dot notation at the end of 
(1571) means the function after the operator and this last term corresponds to 
the discounting. 

C,^i^V2p,Yf,(y)S,-—---/2vA(y)— (56) 



4 i/^ (y) Sf^ + n (^^S. A - .) ^ (/(.» (57) 

If the same procedure of section [2] is followed, equations (l58ll and ( l59il are 
obtained. These equations give the zero and first order terms of the solution, 
which are now denoted by Pq and Pi. Equation fl58|) is equivalent to equation 
f|T6|) and gives the zero order term Pq whose final condition corresponds to 
the vanilla call = 1) or put {ip = —1) option with strike K. When the 
expectation of £2 with respect to the invariant distribution is calculated, the 
Black Scholes operator CBsi^'i) is obtained. This means that the zero order 
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solution is the vanilla Black Sclioles price with volatility equal to (Xj. To 
obtain the first order term Pi, consider equation (160 p which is equivalent to 
equation fl25|) . Now, there is only one term which corresponds to underlying i 
(the other underlying and the correlation term are not present). The function 
(pi is the solution of a Poisson equation similar to the first equation of (j24|) . 

Po{T,S.) = [ifiS,-K)]+ 
£2) Pi = UiC,' (£2 - (£2)) ) Po (59) 



^0 ' (^2 



A))Po = ^0.(l/)5f|^ (60) 



Applying the operator Ci to equation fl60|) . taking expectations with re- 
spect to the limit distribution of Cq and multiplying by y/e yields equation 
( 16T|) . where Ri and Ui are given by equation ( l62i) . See that when the operator 
£1 is applied, the only function which depends on y is (pi{y)- 



CiC^^ (£2 - {£2))) = RA^ i Sf^] + U.Sf^ 



2Ri + Ui) Sf^^ + RiSf^T^ 



dSi \ '^'1. dSf I ' ^^'^^ dSf 

i ds'f ' as? 



(61) 



= C. = -^.„(A^;) (62) 

If equation (159!) is multiplied by y/e, the expectation of £2 with respect 
to the invatiant distribution is replaced in the left hand side of equation (l59i) 
and equation (!6T]) is replaced in the right hand side of equation (!59|) . equation 
( 1^ is obtained. This equation is the same as (5.36) of ^. Identifying terms 
yields V2 and V3 as defined in [2] in terms of Ri and Ui according to equation 

(ED. 

Cbs (ct.) (v^Pi) = {2R, + Ui) + R^S^^ (63) 

\/2 = 2Pi + Ui Vs = Ri (64) 

It is easy to check that the solution of equation (l63l) (the vanilla option 
price) is given by equation fIBSl) (considering that applying Cbs on any of the 
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derivatives with respect to Si gives zero because CbsPq = 0). See that for 
pricing purposes, the market price of risk A{Yt) is set to zero and therefore, 

v/iA (t, S,) = - (T - t) (^{2R, + f/,) + R^S^^^ (65) 

P {t, 5,) = Po (t, S,) + v^Pi (t, S,) (66) 

Equation ( 1671) gives an approximation of the imphed volatihty in terms of 
V2 and V3 where r = St^ fsds and = It^ Qi^^ds. This equation is 
taken from equation (5.55) of section 5.3 of [2] (see that (r — g^) appears now 
instead of r because the process here includes dividends). Setting V2 = 2Ri 
{Ui is equal to zero as the market price of risk is zero) and replacing V3 = Ri 
yields the final calibration equation (168|) . where FiT = Si{t)e^'^~'^^^^'^~''^^ is the 
forward price of underying i at expiration. 
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